Superconducting proximity effect in a diffusive ferromagnet with spin-active interfaces 
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We reconsider the problem of the superconducting proximity effect in a diffusive ferromagnet 
bounded by tunneling interfaces, using spin-dependent boundary conditions. This introduces for 
each interface a phase-shifting conductance G</> which results from the spin dependence of the phase 
shifts acquired by electrons upon scattering on the interface. We show that strongly affects the 
density of states and supercurrents predicted for superconducting/ferromagnetic hybrid circuits. We 
show the relevance of this effect by identifying clear signatures of G<f, in the data of T. Kontos et al 
[Phys. Rev. Lett. 86, 304 (2001), ibid. 89, 137007 (2002)]. 

PACS numbers: 73.23.-b, 74.20.-z, 74.50.+r 
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Superconducting/ferromagnetic (S/F) hybrid struc- 
tures raise the fundamental question of what happens 
when two phases with different broken symmetries inter- 
act. When a F metal with uniform magnetization is con- 
nected to a BCS superconductor, the singlet electronic 
correlations characteristic of the S phase propagate into F 
via Andreev reflections which couple electrons and holes 
with opposite spins and excitation energies. In the dif- 
fusive case, this propagation occurs on a scale limited by 
the ferromagnetic exchange field. The decay of the corre- 
lations in F is accompanied by oscillations of the super- 
conducting order parameter because the exchange field 
induces an energy shift between the correlated electrons 
and holes Q, Q • This has been observed experimentally 
through oscillations of the density of states (DOS) in F 
3, or of the critical current Iq through S/F/S structures 
IHIIQ, with the thickness of F or the temperature. 
Remarkably, the oscillations of Iq have allowed to ob- 
tain 7r-junctions, i.e. Josephson junctions with Iq < 
H, which could find applications in the field of super- 
conducting circuits [2j- 

The interface between a ferromagnet and a non- 
magnetic material can scatter electrons with spin paral- 
lel or antiparallel to the magnetization of the ferromag- 
net with different phase shifts. The spin- dependence of 
the interfacial phase shifts (SDIPS) is a general concept 
in the field of spin-dependent transport. The SDIPS im- 
plies that spins non collinear to the magnetization pre- 
cess during the scattering by the interface. This so-called 
spin mixing is expected to affect drastically the behav- 
ior of F/normal metal when several F electrodes 
with non collinear magnetization are used, as observed 
experimentally by [Tll |. The same phenomenon is pre- 
dicted to occur in F /coulomb blockade island 0, and 
F/Luttinger liquid hybrid circuits. In S/F hybrid 
systems [lj, |lja, |lj| , the SDIPS is even predicted to af- 
fect the system in collinear configurations, due to the 
coupling of electrons and holes with opposite spins by 
the Andreev reflections. However, few experimental sig- 
natures of the SDIPS have been identified up to now (e.g. 
Ref. 0| proposes for the data of 0] an interpretation 



based on the SDIPS). 



In this Letter, we reconsider the problem of the super- 
conducting proximity effect in a diffusive F. Up to now 
the tunnel S/F contacts used to produce this effect were 
described (see e.g. Q) with spin independent boundary 
conditions (BC) derived in 18] for the spin-degenerate 
case. Instead of that, we use spin-dependent BC based on 
Ref. ^(J. These BC introduce a phase shifting conduc- 
tance G,f, which takes into account the SDIPS. We show 
that G<£ strongly affects the phase and the amplitude of 
the oscillations of the DOS or Iq with the thickness of 
F. Our approach thus provides a framework for future 
work on S/F diffusive circuits with tunneling interfaces. 
We show its relevance by a comparison with the data of 
0, 0] which shows that strong experimental manifesta- 
tions of the SDIPS have already been observed through 
the superconducting proximity effect. 



We consider a S/F hybrid circuit with a single F elec- 
trode homogeneously magnetized in direction z . In 
the diffusive limit, the electrons in a superconduct- 
ing or ferromagnetic electrode a can be described with 
quasi-classical and diffusive Green's functions G a in the 
Keldysh<g>Nambu<g>spin space (we use the notations of 
(lif). The BC at a S/F interface can be calculated by 
assuming that the interface potential locally dominates 
the Hamiltonian, i.e. at a short distance it causes only 
ordinary scattering (with no particle- hole mixing). We 
characterize this scattering with transmission and reflec- 
tion amplitudes t„,^ and rf ^ for electrons coming from 
the S(F) side in channel n with a spin a parallel (a =|) or 
antiparallel (a =1) to z . In practice, the planar S/F con- 
tacts used to induce the superconducting proximity effect 
in a diffusive ferromagnet are likely to be in the tunnel 
limit 0,H3j due e.g. to a mismatch of band structure 
between S and F, thus we assume T„ = J2a \^n,a\ 2 ^ 1- 
We also consider that the system is weakly polarized. 
Following HiHHI, the BC at the right hand side F of a 
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FIG. 1: Zero energy density of states at x = d in a S/F/I 
structure, in terms of [(JV(0)/iVo) — l]/7t as a function of 
d/^F, for 7t -C 1 and different values of 7^ (full lines). The 
dotted lines show 4[(N(0)/No) — l]/7t at 1 = d in a semi- 
infinite S/F structure with the same values of y t and 7^,. The 
inset shows the DOS at x = d as a function of 7^ for the 
S/F/I structure. 



S/F interface is 
dG F 



2g F G! 



dx 



G t G s + tG^a z f 3 + ^*D+,G F 
[iG x G s D- + iG^D-Gf, G f ] 



(1) 



with D± = <J z fy,Gs ± Gs<J z T3- Here, a z and T3 are 
Pauli matrices in spin and Nambu space respectively. 
The conductivity of F times area of the junction, noted 
gF, is assumed to be spin independent. The coeffi- 
cient G t = Gq J2 n T n is the tunneling conductance, 
G M r = Gq I]„(|if iT | 2 - I^.J 2 ) is the magnetoresistance 
term which leads to a spin-polarization of the current, 
and G = 2G Q J2 n {Pn ~ 4 [r^/T„]) is the phase-shifting 
conductance, with r,f = im[^ T ^*], p F = Im[r£ T r£j*] 
and Gq — e 2 /h. These three terms already appeared in 
(lif for studying normal electrodes in contact with S and 
F reservoirs (with no proximity effect in F). The extra 

terms in G £ = -Gq£„ t 7? and G x = G eE n T "(P« + 
r^)/4 occur because there are superconducting correla- 
tions at both sides of the interface. Note that G^, G x and 
Gj can be finite only if the phase shifts acquired by the 
electrons upon reflection or transmission at the interface 
are spin-dependent. The exact values of these conduc- 
tance coefficients depend on the microscopic structure 
of the interface. However, we can estimate their rela- 
tive orders of magnitude in a rectangular potential bar- 
rier model, by describing the ferromagnetism of F with 
an exchange field E ex much smaller than the spin av- 
eraged Fermi energy Ep of F. This gives expressions of 
Gmr, G%, G£ and Gj? linear with E ex /E F - The tunnel 
limit can be reached by considering a strong mismatch 
between the Fermi wavevectors in S and F (case 1) or 
a high enough barrier (case 2). In both limits we find 



\Gmr \ , \ G X \ , |Gf| <§; G t , which allows us to neglect these 
terms in the following. In case 1, we find \G<p\ <C Gt 
whereas in case 2 |G<^| can be larger than Gt- Thus we 
will study the consequences of the spin dependent BC for 
an arbitrary value of \G^ \ /Gt- In addition, in case 1 we 
find Gtf, < but in case 2 the sign of G^ depends on the 
details of the barrier, thus we will consider both signs for 
G^. 

In equilibrium, we can use normal and anomalous qua- 
siclassical Matsubara Green's functions parametrized re- 
spectively as cos(A (T ) and sin(A cr ) exp(iip a ) to describe 
the normal excitations and the condensate of pairs (see 
e.g. [13] )• The spatial variations of the superconduct- 
ing correlations in F are described by the Usadel Eqs. 
dQa/dx = and d 2 A a /dx 2 = fc^sgn(w n ) sm(A a )/£ F + 
Q 2 rCOS (A ff )/sin 3 (A CT ), with £p = (hD/E^) 1 / 2 , w„ = 
(2n + 1) ttIcbT. Here, Q a = sm 2 (A rT )d<f rT /dx is the spec- 
tral current (constant with x) and D the diffusion coeffi- 
cient. We introduced k a = (2 (iasgn(uj n ) + \uj n \/ Eex)) 1 ! 2 
for later use Q . Neglecting Gmr , G x and G{ in JQ) yields 

8A„ 

iG^a sm(A a ) + G t [cos(A s ) sin(A CT ) 



gF- 



dx 



sin(A s ) cos(A CT ) cos(yv - tp s )] 



(2) 



gF- 



dx 



sin(A (T ) = G t sin(As) sin^ - ip s ) (3) 



In Eqs. JU and ©, we used rigid BC for S, i.e. A a = 
As = arctan[A/w„], with A the gap of S. 

In the following, we consider the limit of a weak prox- 
imity effect in F, i.e. A CT = 6 a for uo n > and A CT = tt — a 
for w„ < with |6* CT (x)| <C 1. We first study geometries 
with Q a = 0, i.e. no supercurrent flows through the 
device. In this case, the proximity effect in F can be 
probed through measurements of the density of states 
N(e) = N (l - £ a Re[0 2 (:r)]/4) (with u n = -is + 0+ 
and sgn(o;„) = 1). The simplest case of a single S/F 
interface with F at x > yields 

7t sin(A s ) 



aSF 



(x) 



7t |cos(A s )| + «70CTSgn(w„) + k, 



— k - 



(4) 



with 7t(0) = Gt^^F I 9f- In the limit A <C E ex where 
k a = 1 + i(7Sgn(oj n ) , the weak proximity effect hypothesis 
leading to 10} is valid for any values of 7^ and e if 74 <C 
1. Since k a has finite real and imaginary parts, 9^ F (x) 
shows the well-known exponentially damped sinusoidal 
oscillations with d. The remarkable point in is that 
70 shifts these oscillations and modifies their amplitude 
[see Fig. ^ which shows the DOS following from Q]. 
We also study the S/F/I geometry, with F at x € [Q,d] 
and the insulating layer I at x > d, for later comparison 
with the experimental data of |3|. Using for the S/F 
interface and d8 a /dx = for F/I yields: 



3 s a FI {x)=e d a coshKa; 



d) k f 



cosh 



ka ~z 



(5) 
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with 9 c j = 7 t sin(As)/(7t|cos(A s )| + ij,pasgn(uj n ) + 
k a twh{k a d/^F))- In the limit A <C E ex and d > 
the 9~ linearization leading to © is again valid for any 
70 and e if j t <C 1. From Fig. ^ l<j> has qualitatively 
the same effect on 9^ FI (x) as on 9^ F (x) . More quantita- 
tively, for d > £ F one has 9 s / 1 {x = d)/9 s /{x = d) = 2 
|23| and for lower values of d, this ratio depends on d. 

Another way to probe the superconducting proxim- 
ity effect in F is to measure the supercurrent through a 
S/F/S Josephson junction. We consider a junction with 
F at x € [0, d] and a right(left) superconducting reservoir, 
called R(L) at a constant phase (—)(fs/2, A supercurrent 
Is = 7r 9Fk B TJ2nez,a=±iQ^( UJ n)/2e flows through this 
device 0- We focus on the asymmetric limit 7^ <C 7 t L , 
which corresponds to the experiment of H, and assume 
70 = [24]]. We allow L and R to have different su- 
perconducting gaps A R ( L \ so that A CT = Ag L ^ in R(L). 
Solving this problem perturbatively with respect to the 
S/F/I case yields 
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FIG. 2: Critical current Jo of an asymetric S/F/S junction as 
a function oid/^F, calculated from Eq. @ for 7^ <C "ft ^ 1, 

A HR) = A < ^ fesT / A _ 15 ( fuU lmeg ) The 

dashed lines show the large d/^F approximation of Eq. 0. 
The inset is a phase diagram indicating the equilibrium state 
of the junction (0 or n) depending on 70 and d/^F- 



Q CT (w n ) = ^7*sin(A«)sm(<^) 



£f cosh 



k — 



1 -1 



(6) 

where 9^ corresponds to the expression given above with 
As = Ag and 7t(0) = lt(4>)' ^ ne supercurrent has the 
form Is = Iq sm(tps) because most of the phase drop 
occurs at R. In the limit A L = A R = A < E ex , 7 t L < I 
and d/^F > 1, © yields 



el 



= 7r tanh 



2fc B T 



(7) 

with A (7^) = arg[i — (1 +7^)]. It is already known that 
the state of the junction depends on d. Equation (J7J) 
shows that 7^ shifts the oscillations of the Io(d) curve. 
Thus, for a given value of d, the state of the junction can 
be as well as tt, depending on 7^. Fig. [21 shows that 
this effect still occurs when one goes beyond the large 
d/^p approximation. Note that in the limit A <C E ex 
and 7 t L <§; 1 used to obtain Q), it is not possible to find 
a temperature cross-over for the sign of Iq as observed 
in 0, El • However, we expect to find such a temperature 
cross-over with a 0/ir or tt/0 transition, depending on the 
value of r ]i, if the energy dependence of k a is taken into 
account |25| . 

To show the relevance of our approach, we compare our 
predictions with the measurements of Refs. (HTH . We 
first consider the |io| measured in an asymmetric S/F/S 
junction, i.e. Nb/Pdi-^Ni^/Alox/Al/Nb with x ~ 0.1 
and 7^/7^ ~ 10 5 0- We assume that the contacts have 
r„ « 1, which allows to use Eqs. J2EJ- We will use 
the experimentally determined values A Al / Nb = 0.6 meV 
and A Nb = 1.35 meV<C E ex , which implies k a ~ 1 + ia, 
and T = 1.5 K. Samples with different thicknesses d of 



PdNi were measured (see Fig. |3J). Interpreting these 
data requires a careful analysis of the influence of d 
on the different parameters. We have gF = 2e 2 NoDA 
and £p = y/ hD /E ex , with D = vpl/3 and A the con- 
ductors cross-section. Curie temperature measurements 
show that the exchange field E ex increases linearly with 
d 26]. In addition, we first assume that the mean free 
path / is constant with d, as confirmed by resistivity mea- 
surements for d > do = 80 A. This allows to parametrize 
the problem with 7 t L = a° £ F / \/do~d and £f = ^ F \Jdo/d 
where £ F and a° = G L do/gF are constant with d. We 
also assume that G L is proportional to E ex as found 
above in the rectangular barrier model for E ex -C Ep, 
so that we take 7^ = 7° yd/do with 7° constant with d. 
We neglect 7^ due to the existence of a strong insulating 
barrier at R |2J] ■ The absolute amplitude of E ex was not 
determined exactly, so that £ F can be considered as a fit- 
ting parameter as well as a° and 7?. This makes in total 
three fitting parameters but we expect to find for a° a 
value close to the value 0.2 found from minigap measure- 
ments in Nb/Pd [2^]. We have calculated |Jo| by sum- 
ming I© on energy and spin. It is not possible to account 
for the data with 7° = 0. On the contrary, a good agree- 
ment with the experiment is obtained by using q P = 0.4, 



(,° F = 36 A and 



4 



-1.3 (full lines in Fig. 0) 2J1 12 



We have checked that this choice of parameters fulfills 
the hypothesis |^ CT (^)| "C 1 made in our calculations. Re- 
markably, for d ~ do in Fig. the theory for 7° = —1.3 
gives Iq < in agreement with subsequent experiments 
[30j, whereas it gives Io > for 7° = if one keeps the 
same orders of magnitude for a" and £ F . For d < do, I is 
linear with d, which we have taken into account by using 
£f = 70 = 707 a nd 7 t L = a^ F /d, with the same val- 
7° and £ F as previously (dashed lines in Fig. 
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FIG. 3: Critical current measured by [f| through 
Nb/Pdi-zNix/Alox/Al/Nb junctions as a function of the 
thickness d of Pdi-zNLc (symbols). The lines are theoreti- 
cal curves calculated from Eq. @ for d > do (full lines) and 
d < do (dashed lines), with the fitting parameters a t = 0.4, 
£,p — 36 A and the experimentally determined parameters 
A^ 6 = 1.35 meV, A Al/Nb = 0.6 meV, d = 80 A and 
T = 1.5 K. The data are well fitted with 7^ = -1.3. We 
also show the theory for 7° = 0. Inset: DOS measured by 
Q in Nb/Pdi-zNiz/Alox/Al, as a function of d. The full 
and dotted lines show the DOS at x = d calculated from the 
second-order generalization of © (see text), for d > do and 
d < do respectively. We used = 50 A and 7^ = —1.6 or 
7$ = 0, all the other parameters being unchanged. 



|3J. This approach gives a surprising agreement with the 
data, which seems to indicate that the Usadel descrip- 
tion still works for d < do although I is linear with d [31| . 
Kontos et al. have also performed DOS measurements 
in Nb/Pdi-xNiz/Alox/Al 0, prior to the I measure- 
ments |27| . We have assumed again that E ex was linear 
with d in these measurements, to try to interpret the 
N(0) = f(d) curve with the same fitting procedure as for 
io- We have generalized Eq. JSJ) to second order in 9 a 
because the values of d/£p are slightly lower than for the 
io measurements. Again it is impossible to interpret the 
data with 7" = 0. We obtain a satisfactory fit by choos- 
ing £f = 50 A and 7° = —1.6, all the other parameters 
used being the same as in the previous case. Finding a 
£p higher than for the Iq data is in agreement with the 
fact that the samples used for measuring the DOS were 
realized with a lower concentration x of Ni. 

In summary, we have studied the effect of spin- 
dependent boundary conditions on the superconducting 
proximity effect in a diffusive ferromagnet bounded by 
tunneling interfaces. We have shown that the phase- 
shifting conductances G^, describing the spin-activity of 
the interfaces in this context, strongly affect the behavior 
of the system and allow a consistent microscopic expla- 
nation of the DOS and supercurrent data of 0, 0] This 
suggests that such effects will have to be considered in 
any future work on S/F hybrid circuits. In the context 



of spintronics, this approach might also provide a way to 
characterize spin-active interfaces. 
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